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ABSTRACT 


A FORTRAN computer code (SKINTEMP) has been developed to calculate 
transient missile/aircraft aerodynamic heating parameters utilizing basic flight 
parameters such as altitude, Mach number, and angle of attack. The insulated skin 
temperature of a vehicle surface on either the fuselage (axisymmetric body) or 
wing (two-dimensional body) is computed from a basic heat balance relationship 
throughout the entire spectrum (subsonic, transonic, supersonic, hypersonic) of 
flight. This calculation method employs a simple finite difference procedure which 
considers radiation, forced convection, and non-reactive chemistry. Surface 
pressure estimates are based on a modified Newtonian flow model. Eckert’s 
reference temperature method is used as the forced convection heat transfer model. 
SKINTEMP predictions are compared with a limited number of test cases. 
SKINTEMP was developed as a tool to enhance the conceptual design process of 
high speed missiles and aircraft. Recommendations are made for possible future 
development of SKINTEMP to further support the design process. 




n \fu>it> 

CMIW 

THESIS DISCLAIMER 

The reader is cautioned that the computer program SKINTEMP 
has not been exercised for all cases of interest. While every 
effort has been made within the time available to ensure the 
program is free of computational and logic errors, it cannot be 
considered validated. Any application of this program without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. GENERAL 

A FORTRAN computer program has been created to determine 
transient missile and aircraft mission aerodynamic heating 
parameters such as heat transfer coefficients, insulated skin 
temperatures and heat flow rates at some point on the surface 
a wing or a body of revolution. Required data supplied from 
the user include material properties and trajectory data such 
as Mach number, angle of attack, altitude, and time of flight. 
This program calculates aerodynamic heating parameters for 
subsonic, transonic, supersonic and hypersonic (no reacting 
chemistry) flight regimes. 

B. BASIC THEORY 

The governing equation which serves as the basis for the 
computer program incorporates both radiative and forced 
convective heat transfer processes. The thermal energy 
balance equation for this insulated skin heat transfer case is 
given by 

= hA(T aw - T e ) - oeA(T/ - T r 4 ) . (1) 

The equation written in this form assumes that the adiabatic 
wall temperature, T„„ is greater than the insulated skin 
temperature, T s , which is greater than the radiation reference 




temperature, T r . The heat transfer coefficient for convection 
is obtained from the Nusselt number which is a function of the 
Reynolds and Prandtl numbers. These dimensionless numbers are 
defined as 

(a) = 2 ™ (jb) N = hx ( 2 ) 
k " k 

Once the complex heat transfer coefficient calculation is 
performed and the other significant heat transfer parameters 
determined, a finite difference numerical method can be 
applied, yielding a skin temperature profile for the specified 
point of interest. 

C. SCOPE 

Chapter II details the theory used to develop the SKINTEMP 
program. 

Chapter III describes the actual computer program, its 
required entry data, expected output data and general program 
limitations. 

Chapter IV presents test case results with the appropriate 
comparisons, which include intermediate heat transfer 
calculations. One test case is obtained from a heat transfer 
textbook and compared to corresponding SKINTEMP results. A 
less descriptive NASA test case is also compared with SKINTEMP 
results. 

Chapter V provides recommendations for future code 


development. 



II. 


TRANSIENT HEAT FLOW THEORY 


The fundamental concepts used to build the theoretical 
model for SKINTEMP involve two basic modes of heat transfer, 
convection and radiation. Forced convection is assumed because 
the vehicle is propelled through the air by the release of 
chemical energy. During supersonic and hypersonic flight, 
the local stagnation pressure at the edge of the boundary 
layer is assumed to be defined by the stagnation pressure 
behind a normal shock wave. Figure 1 illustrates this 
concept. 



Figure 1. Stagnation pressure at boundary edge 







Boundary layer flow is assumed to be turbulent due to the 
magnitude of the Reynolds number which is generally greater 
than 500,000. The gas dynamic relations are based on inviscid 
flow, no reacting chemistry assumptions [Ref. 1]. 

The primary objective of the computer code is to rapidly 
compute the transient insulated skin temperatures along the 
trajectory of the aircraft at a given point somewhere on the 
surface. To accomplish this, the thermal energy balance 
equation 

mC D 4^ = hA{T - T ) - oeA(T. 4 - T r 4 ) (1) 

* o t 

is solved for T„ using a finite difference method [Ref. 2], 
However, before this can be accomplished, the adiabatic wall 
temperature and the heat transfer coefficient must be 
determined. The required data for the skin temperature 
calculations include trajectory parameters (time, Mach 
number, angle of attack, altitude), material properties 
(density, material specific heat, emissivity), the location 
and skin thickness at the point of interest, and a radiation 
reference temperature. This radiation temperature is the 
temperature of a distant surface seen by the insulated skin 
element, which for this program, is either space or the earth. 



THE ADIABATIC WALL TEMPERATURE 


The first parameter needed to determine skin temperature 
is the adiabatic wall temperature. The adiabatic wall 
temperature is given by 

T„ = T l [l + . (3) 
The local Mach number needed for the solution. 



is dependent on the local pressure and the local stagnation 
pressure. The speed of the aircraft, whether subsonic or 
supersonic will determine the local stagnation pressure. If 
the trajectory profile calls for subsonic flight, then the 
local stagnation pressure will be 

P tL = ?l h + ^Mf] Y/T_1 . (5) 

If the Mach number is in the supersonic range, the local 
stagnation pressure is then 




P tL 


Ml 


. + 


y+lj 


(6) 
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Freestream pressure, temperature, and density are calculated 
using standard atmosphere property value approximations found 
in Reference 3. 


Local surface pressure is determined using modified 
Newtonian theory, [Ref. 4] where 

P L = CpnaxCOsZQq + P m (8) 

p„ r - p 

- — (9) 


(10) 


g 

g = l/2p„vf. 

Modified Newtonian flow theory has been shown to be applicable 
for the prediction of local surface static pressures over all 
surfaces experiencing non-separated flow. [Ref. 5,6] The 
accuracy of Modified Newtonian flow theory is quite reliable 
at hypersonic speeds. Its use at subsonic and low supersonic 
speeds is justified by the continuity and simplicity it 
provides in the lower Mach number regions where aerodynamic 
heating rates are so low as to be negligible. 

The remaining parameter needed to determine the adiabatic 
wall temperature is the recovery factor (turbulent flow is 
assumed), given by 

r = N& (H) 


where the Prandtl number is 


N pr 


C p Vi 

k 


( 12 ) 


The specific heat of air used in the Prandtl number is given 
by [Ref. 2] 

C p = .219 + .342£’-4T t - .293£-8Tf. (13) 

6 



Thermal conductivity [Ref. 7] and dynamic viscosity [Ref. 2] 
are also temperature dependent and are given respectively by 



2.21E-S Tl /2 
T L +198.7 


B. THE HEAT TRANSFER COEFFICIENT 

In order to calculate the heat transfer coefficient,one 
uses the concept of the reference temperature. The reference 
temperature is a temperature whose magnitude is typically 
between that of the wall temperature and freestream 
temperature and in low speed flow, the reference temperature 
does become the average of the two [Ref. 2]. This reference 
temperature more accurately models the temperature 
distribution that occurs within the boundary layer, which has 
a significant impact on the skin temperature [Ref. 8] . The 
transport properties should be evaluated at the reference 
temperature [Ref. 9] which is defined by 

T* = T l + .5 {T s - T l ) + . 22 ( T aw - T l ) . (16) 

Values for the thermal conductivity, dynamic viscosity, and 
specific heat used in the Nusselt number are based on this 
reference temperature. 

The Nusselt number can be defined as 

Nm = 

7 


h X 





Based on experimental data, the Nusselt number can also be 
written as a function of Reynolds and Prandtl number, [Ref. 2] 
Nm = .0292 (N^) e (Nn) 1/3 (18) 

where the boundary layer is assumed to be turbulent and the 
Reynolds number is defined as 


Nr E 


p'V L X 

n* 


(19) 


With these significant variables defined, the thermal 
energy balance equation can be solved using a simplified 
finite difference method. This method calculates the 
temperature value at a future time based on a calculated time 
increment and the current temperature value. 


C. THE SKIN TEMPERATURE 

Solution of the thermal energy balance equation for the 
incremental skin temperature yields 

T s i ,i = [ hA(T aw - T a ) - aeA(T s 4 - T r 4 )]-^. + T sl . (20) 


For stability purposes, the time interval associated with this 
finite difference method is given by 




At 


mC p 

*T+h~ r 


h r = Aeo [rj + r|r r + tIt s + rj] . 


(21) 


( 22 ) 


Each progressive skin temperature is dependent on the 


previous temperature, heat transfer coefficient and adiabatic 



st meet the stability requirement of the finite difi 
thod. This stability requirement, which preven' 
'lution from diverging, is given by 


To meet this requiremei 


! interval was chosen ; 


D. THE HEAT FLOW RATE 

The heat flow rate is dependent on the location of the 
stagnation heating point. In the stagnation region for the 
axisymmetrical and two dimensional platforms, the maximum 
stagnation point is 70 degrees from the centerline [Ref. 10] 
and the stagnation heating rate is given by [Ref. 11] 

a = 17600 r~p~ r u t 15 [ Hsi ~ h ” i (25) 

sta9 Po KJ [h s1 - H w545 \ 

Away from the stagnation point but prior to regions that can 
be treated as flat plates, the heat flow rate [Ref. 10] is 


as where flat pla- 


g heac = hA(T at 


(27) 





III. THE SKINTEMP COMPUTER PROGRAM 


A. THE COMPUTER CODE STRUCTURE 

The flow logic used for the SKINTEMP program is 
illustrated in Figure 2. The program is divided into four 
major sections, data entry, atmospheric property generation, 
skin temperature calculations and data output. 

The first module contains the data entry loop. The data 
can be entered via two methods, keyboard entry or file entry. 
If the user chooses to utilize keyboard entry, the program 
will prompt the user for the specific data. Each entry should 
be followed by a keyboard hard return. The file entry method 
requires the user to create an external file using any 
standard editor containing all the required information. For 
file entry, the correct syntax to enter at the prompt is 
SKINTEMP. in. The .in refers to the entry file; any input file 
used must have the .in included in the file name. This 
command allows data to be read automatically by the input 
subroutine. 

The important note to remember is that for file entry, 
each piece of fixed point data must be placed on a separate 
line and the entry data must correlate with the data requested 
by the program. If the program terminates with an error 
message, this would be the first place to check for mistakes. 
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The second section of the code generates standard 
atmospheric values that are needed for the freestream 
pressure, temperature, and density calculations. The 
governing equations used to model atmospheric properties are 
divided into two regions, gradient and isothermal [Ref. 3], 
The gradient relationships are 


p T -£ 9/LR *i] 

Po (28) 

_P = (29) 

P n T a 

(30) 

T = L(h-h ) + T a . 


There are several isothermal regions in the atmosphere. 
While the temperatures remain constant in each of these 
isothermal regions, the temperatures of the gradient regions 
will vary. [Ref. 3] The pressure and density relationships 
for isothermal regions are given by 


_P_ 

P Q 


e -3p-h 0 ]/RT 


(31) 


_e_ = e -S[A-Ao]/ KT (32) 

Po 


With these relationships, the atmosphere is modeled to 260,000 
feet with isothermal layers from approximately 36,000 to 
82,000 feet (temperature = 392.4 °R, density = .00071 slug/ft 3 , 
pressure = 476 lb/ft 2 ) and 154,000 to 174,000 feet 

(temperature = 509 “R, density = .00000303 slug/ft 3 , pressure 
= 2.6419 lb/ft 2 ) . 
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The third major section of SKINTEMP is the calculation of 
the skin temperature. This section is a large counted loop 
because each successive skin temperature is dependent on the 
previous temperature calculated. This loop is partitioned 
into four sections. The first section calculates adiabatic 
wall temperature and the second determines the heat transfer 
coefficient. These two calculations lead to the skin 
temperature solution, which is in the third section. The 
fourth and final section calculates the heat flow rate. [Ref. 
Chapter II] Appendix A contains the actual computer program. 
The initial documentation of the program contains a list of 
the variables, their definitions, and the English units used. 

It is important that the user understand the differences 
between the time values used throughout SKINTEMP. The user 
should input time intervals that best linearize the 
trajectory. In other words, whenever the Mach number, 

altitude or angle of attack is changing, this time and 
associated flight parameters should be a data entry point. 
These time intervals can be the print out time intervals but 
SKINTEMP also allows the user to print the results at 
different time intervals if desired. These input and output 
time intervals are independent of the time interval chosen for 
the heat balance calculations. Figure 3 demonstrates these 
differences and should be used as a reference when developing 
a data entry strategy for SKINTEMP. 


12 



TAW IS COMPUTED USINI 
P tL , V_, q, Q l , C pmax , P L , 
T tL , k, C p , |X, N pr , r. 


T r ,p z ,Vi z ,V L ,N 
C p , k r , Npn, N t 



COMPUTED AFTER CALCULATING 


LOW RATE COMPUTED BASED ON STAGNATION POINT 


£ 6 t < INPUT TIME INTERVAL 


£ 5 t * INPUT TIME INTERVAL 


Figure 2. SKINTEMP Flowchart 
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Figure 3. Time Interval Comparist 





REQUIRED DATA INPUT 


SKINTEMP allows the user to enter data through 
the keyboard or through an external file. In either case, the 
following data is required to successfully run the program: 


• The shape of the body, which must be either two 
dimensional or axisymmetric 

• If a two dimensional wing shape is selected, the wing 
sweep (degrees) 

• The launch surface skin temperature of the vehicle (°F) 


• For each time interval, the Mach number, the angle of 
attack (degrees), and the altitude (feet) 


WARNING: Entering a Mach number of 
zero will cause the 
program to abort. See 
Section D. below. 


• The distance from the nose or leading edge to the point of 
interest (feet) 

• The thickness of the surface at the point of interest 
(inches) 

• The nose radius (feet) 

• The material density (pound per cubic foot) 

• The material emissivity (dimensionless) 

• The material specific heat (BTU per pound-°F) 

• The orientation of the surface. Does the surface face the 
earth or space for the majority of the flight? This 
provides the radiation reference temperature of -460°F for 
space and 60 °F for the earth. (°F) 

• The angle from the unit surface normal of the skin element 
to the freestream (degrees) 

• The output time interval (seconds) 
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C. SKINTEMP OUTPUT 

The output of the computer program is in tabular form. 
The user can choose any time interval for printing independent 
of the data time entries. See Appendix B for an sample 
output. SKINTEMP will also produce an external UNIX file 
named Fort.3 which produces a table of values for the time and 
skin temperature that can be used in typical engineering 
graphics/plotting programs. 

The program is clearly documented and the print subroutine 
is the last section in the program. The output format is as 
follows: 

TIME ADIABATIC WALL HEAT TRANSFER SKIN HEAT FLOW 

TEMPERATURE COEFFICIENT TEMPERATURE RATE 

(sec) (“Rankine) (BTU/sec-ft 2 -°F) (°Rankine) (BTU/sec) 

D. SKINTEMP LIMITATIONS 

The purpose of this code is to provide a design tool for 
conceptual missile and aircraft design projects. At the 
present time, SKINTEMP can only calculate the transient heat 
parameters for bodies that are either two dimensional or 
axisymmetric body shaped. See Chapter V for future 
recommendations for code expansion. 

Another limitation and perhaps the most important for user 
information is that at no time can the Mach number entered be 
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equal to zero. The reason for this is that in the calculation 
for C pnux , dynamic pressure is the divisor. C pmal and dynamic 
pressure are given by 

c p,max = ( p tL ~ p J/q (9) where q = l/2pV 2 = pm 2 . (33) 

A Mach number of zero will cause C pmax to be undefined. If 
this occurs, the program will terminate without results. If 
a Mach number of zero is the actual value achieved in the 
trajectory, a small value should be entered, i.e., .001 Mach. 

The last limitation is the array size of the input data. 
At the present time, the array size of data points entered is 
100. If more than 100 data points are desired by the user, 
the array size should be altered in the body of SKINTEMP. The 
array dimension statement is the first to occur after the 
initial documentation and the maximum limit is dictated by the 
version of FORTRAN being used. 
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IV. TEST CASE RESULTS 


A. TRUITT TEST CASE 

The test case used to validate SKINTEMP is based on an 
example skin temperature problem from Robert W. Truitt's text 
Fundamentals of Aerodynamic Heating, [Ref 12] . The example 
missile profiled is axisymmetric in shape and 40 feet in 
length. The skin material used is Inconel steel and the skin 
thickness of the material at the point of interest, 20 feet 
from the nose, is .00504 inches. The flight trajectory data 
are given in Figure 4. 

Trajectory Mach number 



Trajectory Altitude 
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Truitt assumes turbulent flow but does not use specific 
local or reference state conditions. Instead, he assumes 
freestream conditions throughout the calculations. This 
approach is responsible for at least some of the noticeable 
deviation in the comparative results [Ref. Section 2], 

A Runge - Kutta numerical integration method for the skin 
temperature calculation is employed by Truitt rather than the 
finite difference method used in SKINTEMP. This difference in 
approach was one of the primary reasons this test case was 
used. Similar results from two different equivalent numerical 
methods tend to confirm results. 

1. Data Comparisons 

Throughout Truitt's skin temperature calculation 
example several intermediate heat transfer and fluid flow 
parameter values are defined. These computed values include 
freestream static pressure and temperature, Reynolds number 
and the computed skin temperature (based on freestream 
properties). The comparisons of these parameters are listed 
in the following graphs and tables to validate SKINTEMP and 
illustrate some of the theoretical differences between 


methods. 




a. Freestrean Pressure Comparison 

Figure 5 graphically displays the difference in 
freestream static pressures while Table 1 details the exact 
values computed. Truitt's values for freestream static 
atmospheric density are not as close to the NACA standard 
atmosphere [Ref. 13] as the SKINTEMP code calculations and the 
slight divergence in values contributes to the difference in 
the predicted skin temperatures. 


FREESTREAM PRESSURE COMPARISON 



Tims (seconds) 

Figure 5. Freestream static pressure comparison 
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TABLE 1. FREESTREAM PRESSURE 
TRAJECTORY 


TIMETRUITT 
(SECONDS) (LB/FT 3 ) 


1 _ .00237 

~2 .00234 


_5_ .00211 

_6_ .00202 

_7_ .00193 

_8_ .00184 

_9_ .00177 

10 _ . 00170 



| 30 | .00001 


FOR TRUITT'S 


SKINTEMP NACA 

(LB/FT 3 ) (LB/FT 3 ) 


.00238 ,00238 

.00234 .00234 



.00228 _ .00228 

■ 00225 _ .00224 

.00219 _ ,00218 

.00206 _ .00205 

.00178 _ .00178 

.00176 _ .00176 



.00001 | .00001 


I 

















a. Temperature Comparison For Truitt's Trajectory 

Figure 6 and Table 2 compare Truitt's freestream 
temperature to SKINTEMP's freestream and reference static 
temperature. The freestream static temperature calculations 
compare favorably. Again, the divergence in values can be 
partially attributed to differences in atmospheric modelling. 

Local temperatures are compared with the 
freestream temperatures to illustrate the difference of using 
local atmospheric properties instead of freestream properties. 
SKINTEMP uses local conditions for the heat transfer 
calculations while Truitt utilizes the freestream conditions. 
Temperatures near the skin of a hypersonic vehicle will be 
much higher than the surrounding freestream air due to 
boundary layer interactions (i.e. skin friction). Thus, 
calculations using freestream properties are expected to be 
more approximate than those using the loaal conditions. 



Figure 6. Temperature Comparison 
22 





TABLE 2. 


TEMPERATURE COMPARISON FOR TRUITT'S TRAJECTORY. 







































b. Reynolds Number Comparison For Truitt's 

Trajectory 

The difference in theoretical approaches is very 
apparent when comparing the Reynolds numbers in Figure 7 and 
Table 3. Truitt's Reynolds number calculations are very close 
to SKINTEMP's initial calculations, but again, the divergence 
in freestream atmospheric properties clearly has an impact on 
the freestream Reynolds number. The only significant 
consideration is that this difference becomes larger as the 
program builds on data from the previous time step. 

Figure 7 clearly illustrates that the reference 
state Reynolds number is significantly smaller than the 
freestream value. Recall that the reference state utilizes 
the local flow field values at the edge of the boundary layer 
with the transport properties evaluated at the reference 
temperature, and that the Reynolds number is given as 


where 


Nk 


p'VX 


(18) 


p = P/RT 
V = M^fRT 
and 


2.2TE-8 Tl /Z 



(34) 


(35) 


(15) 
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The Reynolds number is inversely proportional to the 
temperature. Use of the local values for temperature, which 
are significantly larger than the freestream values of 
temperature, will cause the Reynolds number to decrease. 

Another parameter that will contribute to a 
difference in Reynolds number calculations is the relationship 
used to determine viscosity, equation 35. Truitt did not 
utilize the same relationship and even though the values are 
of the same magnitude, the difference only adds to the 
deviation in results. 


FREESTREAM AND REFERENCE REYNOLDS NUMBER COMPARISON 



Time (seconds) 

Figure 7. Reynolds Number Comparison 
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TABLE 3. REYNOLDS NUMBER COMPARISON FOR TRUITT'S TRAJECTORY 



























Temperature (Rankine) 


c. Skin Temperature Comparison 

Figure 8 and Table 4 compare the final skin 
temperature calculations between Truitt and SKINTEMP. 


TRUITT SKIN TEMPERATURE PREDICTION COMPARISON 


1500 



20 25 30 


Figure 8. Skin Temperature Comparis 
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SKIN TEMPERATURE COMPARISON FOR TRUITT' 


TRAJECTORY 

































1 . 


Results F] 


Truitt's Trajectory 


The graphs and tables indicate comparable results in 
the two different methods of calculating skin temperature. 
There are numerous reasons for the slight differences in the 
results. The discrepancy in standard atmosphere values leads 
to significant differences in the aerodynamic parameters, 
which cause marked deviations in the heat transfer 
coefficients, and the adiabatic wall and skin temperatures. 

Another major difference is the use of the reference 
state in SKINTEMP and the freestream state in Truitt's test 
case. Using the reference condition instead of the freestream 
state is considered to be the more correct model of what is 
actually occurring aerodynamically. One aspect of this 
difference is manifested in the Reynolds number calculation. 
Despite the significant difference in Reynolds numbers used to 
compute the Nusselt number, the results are nearly the same. 
Recall that to compute the Nusselt number, the Reynolds number 
is raised to the 0.8 power. Large differences in Reynolds 
number between the two methods are minimized in this 
calculation, allowing the skin temperatures to more closely 
correspond with one another. 

Truitt also uses varying skin material specific heats 
while SKINTEMP assumes an average constant specific heat. 
This does not cause a noticeable change in skin temperature 
calculations but does account for a few degrees of difference. 
The choice of time constants is also critical for stability 
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purposes. Truitt uses a time step of two seconds while the 
average time constant step generated by SKINTEMP is 0.3 
seconds. As indicated by Truitt, the smaller the time 
constant, the more accurate the results. Truitt used the 
larger time constant for ease in calculations because they 
were all hand computed. With all other variables being equal, 
the calculation of skin temperatures with the more simplified 
finite difference method versus the Runge-Kutta method is 
expected to give the same degree of accuracy, provided the 
time step is small. 

While the Runge-Kutta integration method is an 
accurate numerical method, Truitt did not have the same 
computer capabilities that are in use today and the current 
standard atmosphere tables are considered to be more accurate 
than those Truitt used. Even though the SKINTEMP results are 
on average nine percent greater than Truittt's results, the 
skin temperature curves are similar with SKINTEMP producing 
the more conservative worst case results. 

A. NASA TEST CASE 

NASA AMES provided an additional test case for SKINTEMP 
verification. The trajectory is typical of a Single-Stage - 
To-Orbit (SSTO) rocket powered vehicle and the point of 
interest is 75 feet aft of the nose. The thickness of this 
point of interest is 0.035 inches and the material at this 
point is Tailorable Advanced Blanket Insulation (TABI). 
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Temperature (Ran 


The] 


re is a body angle of incidence of 13 degrees in addition 
to the varying angle of attack experienced through flight. 
Figure 9 graphically depicts the results and is followed by 
the detailed listing of skin temperatures. 


NASA SKIN TEMPERATURE PREDICTION COMPARISON 



Figure 9. NASA Skin Tempe 


Comparisc 






TABLE 5. SKIN TEMPERATURE COMPARISON FOR NASA AMES TEST CASE 


TIME 

(SECONDS) 

NASA 

(“RANKINE) 

SKINTEMP 

(°RANKINE) 

1.3 

511 

511 

2.61 

571 

570 

3.48 

590 

630 

4.69 

738 

870 

5.78 

1313 

1104 

6.82 

1570 

1449 

7.70 

1821 

1756 

8.71 

2094 

1979 

9.81 

2318 

2262 

10.76 

2463 

2424 

11.79 

2637 

2550 

12.85 

2799 

2645 

13.72 

2912 

2733 

14.72 

3011 

2806 

15.75 

3127 

2864 

16.75 

3200 

2911 

17.75 

3265 

2974 

18.75 

3308 

3019 

19.70 

3267 

3084 

20.78 

3008 

2606 

21.38 

3231 

2133 








1. Results From NASA Trajectory 

Again, the results from the two skin temperature 
calculation methods compare favorably. The NASA test case was 
also based on freestream conditions with surface pressures 
given by tangent wedge/cone theory. While the use of 
freestream conditions is not theoretically correct, there has 
been significant experimental data generated to justify this 
assumption [Ref. 14] . However, these assumptions are best 
used for upper surfaces. SKINTEMP is not limited to this 
assumption and can correctly calculate skin temperatures for 
upper and lower surfaces. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The SKINTEMP computer program was developed as a tool for 
use in the conceptual design of high speed missiles and 
aircraft. It can be used over a wide range of trajectories, 
profiles and configurations and rapidly yields results that 
enable the user to easily investigate changes in material 
characteristics and flight parameters without having to employ 
more complex and cumbersome calculations. As with any 
computer program, the results can only be as good as the data 
entered. This is especially true of SKINTEMP because valid 
results can only be expected if the trajectory is accurately 
modeled. 

Validation of SKINTEMP was achieved using a textbook skin 
temperature calculation example and a more current NASA test 
case. SKINTEMP compared with the NASA and Truitt test cases 
suggest that the SKINTEMP results are reasonable. While each 
of the methods has its own merits, the fundamental theoretical 
base, direct approach and simplicity of SKINTEMP make it ideal 
as an aerodynamic and astrodynamic design tool. 


34 



RECOMMENDATIONS 


SKINTEMP currently calculates the insulated skin 
temperature and heat flow rate at only one surface location. 
The computer code could be expanded to perform a thermal nodal 
analysis that could predict skin temperatures along an 
extended surface. Moreover, temperatures of materials within 
the air/spacecraft affected by the skin temperature could be 
calculated. This would be very useful in the design process 
to initially identify hot spots which, especially in 
hypersonic flight, could be crucial to the overall design. 
SKINTEMP was written strictly for the turbulent flow case 
because this generally produces worst case skin temperatures. 
The code could easily be modified to model laminar flow if 
turbulent flow was not the required assumption. 

Because the skin temperature problem is quite detailed and 
the inputs are extensive, complete trajectory and 
corresponding surface skin temperatures profiles are difficult 
to find in the literature. Additional test cases could be 
developed for SKINTEMP for further validation as well as any 
experimental data that is available. 

The important objective to remember in any significant 
code modification of SKINTEMP is that the program needs to run 
quickly and allow the user to modify any design parameters. 
If the program becomes too complex and time consuming, then 
the usefulness of SKINTEMP is degraded. 
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APPENDIX 


SKINTEMP COMPUTER CODE 
















1/ ft'3) 

lling (lb£/ft*2) 
lng (slug/ ft“3) 
”2) 

f/ff2) 

(Rankine) 









>£ attack (radians) 

1 enthalpy (BTll/lbm) 

»lpy (BTU/lbm) 

:halpy at 540 II (BTU/lbm) 



























T * EXP (-((g / (R * T(HGT))) * 

4) - 154200.))) 

nsT * P(HGT) / PressT 


* (TABLE(HGT,4) - 82021.) 
HGT)/389.97)**(-g/(L2 * R )) 

9 * ((T(HGT)/389.97)** 

L2 * R))+ 1.0)) 


f (-((9 / (R * T(HGT))) * ( 
) - 36089.))) 

?(HOT) / PressT 






sumtime * del tam/( timed! f ) + TABLE(n, 2 ) 


a tmpC - P(n) + sumtime * (p(n + l) - P(n))/timedlf 

atmdens - Dens(n) + sumtlme * (Dens(n+1) - Dens(n)) 
/timed!f 

atmT - T(n) + sumtime*(T(n+1)-T(n))/timedi£ 


if (BODY .EQ. 2) then 

Ma - Ma * cosfsweep * pi/180.) 
end if 

end if 

if (sumtime .EQ. 0 ) then 

Ha - TABLE( n , 2 ) 
atmpc - P(n) 
atmdens - Dens(n) 
atmT - T(n) 
if (BODY .EQ. 2) then 

Ha - Ma * cos(sweep * pi/180.) 
end if 


end if 

60 if (Ha .LE. 1.0) then 

Ptl - atmpc * (1.0 + ((gamma-1.0)/2.0) * Ma*+2) 
c (gamma/(gamma-1.0)) 


else 


Pt - atmpc * (1.0 + ((gamma-1.0)/2.0) * Ma**2)** 

(gamma/(gamma-1.0)) 

Tt - atmT * (1.0 + ((gamma-1.0J/2.0) * Ma**2) 

Ptl - Pt * (((gamma+1.0)/2 * Ma**2)/(1 + (gamma 
-1.0)/2 * Ma**2))**( gamma/(gamma-1.0) ) * 

(2.0*gamma/(gamma+1.0) * Ma**2 -(gamma-1.0)/ 
(gamma+1.0))**(l/(1-gamma)) 




Adiabatic wall tempecatuce calculations 
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q - .5 * atmdens * V**2 

theta - (ANGLE + TABLE( n, 3 ) ) » pi/180. 

CpMax - (Ptl - atmprj/q 

Pi - q*CpMax * (cos(theta))* *2 + atmpr 

Ml - (((Ptl/Pl)**((gamma-1.0)/gamma)-l.0) * (2.0/ 
(gararoa-1.0)))**. 5 

Ttl - atmT * (1.0 + (gamma-1.0)/2.0 * Ma‘*2) 

Tl - Ttl/(1.0 + (gamma-1.0)/2.0 * Ml**2) 

k - (.0141 * (717.0/(Tl + 225.)) * (Tl/492.0) 
**1.5)/3600.0 

Cp - .219 + .342E-4 * Tl - .293E-B * Tl**2 
Cl - 2.27E-8 
C2 - 198.7 

u - Cl * Tl* * 1.5/(Tl + C2) 

Pt - Cp * u * g /k 
re - Pr**(l.0/3.0) 

TAW - Tl * (1.0+re * (gamma-1.01/2.0 * Ml**2) 


oefficient calculations 


Tstar - Tl + .5 * (Tn-Tl) + .22 * (TAW-Tl) 

denstar - P1/(R * Tstar) 

ustar - (Cl * Tstar**l,5)/(Tstar + C2) 

al - (gamma * R * Tl)**.5 

vl - al * Ml 

Restar - denstar * VI * X/ustar 

Cpstar - .219 + .342E-4 * Tstar - .293E-8 * Tstar 
**2 

kstar - (.0141 * (717.0/(Tstar + 225.)) * (Tstar/ 
492.0)**1.5)/3600.0 

tar * g * ustar/kstar 
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NUstar - .0292 * (Restar**.8) * (Prstar 
“(1.0/3.0) ) 


NUstar 


.185 * (Restar) * (Prstar)** 
(i.0/3.0)/((loglO(Restar))** 
2.584) 


end if 


h - NUstar * kstar/X 



hws - Cp * Tn 

hw540 - Cp * 540. 

stagpt - cos(70 * (pi/180.)) 

if (ANGLER .LE. stagpt) then 

qstag - 17600.0/SQRT(RAD) * SQRT(atmdens/Densin) * (V/ 
26000.0)**3.15 * ((hsl-hws)/(hsl-hw540)) 
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